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Abstract 

We deepen the study of the so-called "silver channel" v e — > v T [I] and of its 
relevance to solve some of the ambiguities that can arise in the simultaneous 
measurement of (#13, <5) at the Neutrino Factory by presenting in full detail 
the characteristics of the considered OPERA-like detector and the experimental 
treatment of the different backgrounds and signals. Furthermore, we perform 
a detailed study of the systematic errors associated both with the OPERA-like 
and the magnetized-iron detectors and their effects on the sensitivity. Finally, 
we also apply a refined statistical analysis of the simulated events based on the 
frequentist approach. 



1 Introduction 



The hypothesis of neutrino oscillations 2\ is at present strongly supported by atmo- 
spheric, solar, accelerator |3] and reactor neutrino data. If we do not consider the 
claimed evidence for oscillations at the LSND experiment [5], that must be confirmed or 
excluded by the ongoing MiniBooNE experiment jB], oscillations in the leptonic sector 
can be easily accommodated in the three family Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) mixing matrix Upmns' 
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with the short-form notation s y - = sin 9y,Cij = cos^y. Further Majorana phases 
have not been introduced, since oscillation experiments are only sensitive to the two 
neutrino mass differences Am^ , Am| 3 and to the four parameters in the mixing matrix 
of eq. (Q): three angles and the Dirac CP- violating phase, 5. 

In particular, data from atmospheric neutrinos and from K2K are interpreted as 
oscillations of into v T with a mass gap that we denote by Am| 3 . The corresponding 
mixing angle is close to maximal, sin 2 26* 2 3 > 0.9, and jAm^l is in the range 1.9 
to 3.5 x 10~ 3 eV 2 j7j. On the other hand, the longstanding solar neutrino problem 
has finally been solved by the combination of the SNO data [8 : and recent KamLand 
results jl], isolating the LMA-MSW region |9 as the only viable solution of the solar 
neutrino deficit with v e oscillations into active (z/^, z/ r ) neutrino states. Two allowed 
regions are identified inside the LMA-MSW solution |10j . corresponding to Anig = 
Am 2 2 ~ 7x 10~ 5 eV 2 (LMA-I) and Am 2 2 ~ 1.4x 10" 4 eV 2 (LMA-II). For both solutions 
the corresponding mixing angle (#12) is large (albeit maximal mixing is excluded at 
3a). Finally, a comprehensive three-family analysis (including the negative CHOOZ 
results [TT]) put a bound on #13, sin 2 6*13 < 0.02. 

The planned long baseline experiments ^21-^3 will improve the measurement of 
Am 2 tm (Am 2 tm ~ Am| 3 ~ Amf 3 ) and of 6*23 and measure or increase the bound on 
#13 [H3 H"7j (see also ^H])- This new generation of experiments, however, is only the 
first step toward the ambitious goal of precision measurements of the whole three- 
neutrino mixing parameter space, including the leptonic CP-violating phase 5. This 
long-lasting experimental program consists of the development of some "superbeam" 
facilities (whose combination can strongly improve our knowledge on # 13 , see |19j ) 
and, eventually, of a "Neutrino Factory" (high-energy muons decaying in the straight 
section of a storage ring, thus producing a very pure and intense two-flavor neutrino 
beam [2D1IS])- One of the main goals of the Neutrino Factory program (see for example 
[2212111 and refs. therein) would be the discovery of leptonic CP violation and, possibly, 
its study |231-|27|. 

The most sensitive method to study this topic is the measure of the transition 
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probability u e (u e ) — > v^iv^. This is what is called the "golden measurement at the 
Neutrino Factory" |2Z1- At the Neutrino Factory an energetic electron neutrino beam 
is produced with no contamination from muon neutrinos with the same helicity (only 
muon neutrinos of opposite helicity are present in the beam, contrary to the case of 
conventional beams from pion decay). Therefore, the transition of interest can be easily 
measured by searching for wrong-sign muons, i.e. muons with charge opposite to that of 
the parent muons in the storage ring, provided the considered detector has a good muon 
charge identification capability. However, the determination of (613,8) at the Neutrino 
Factory is not at all free of ambiguities, much as it was the case for the different 
solutions to the solar neutrino problem before KamLand. In [2E] it was shown that, 
for a given physical input parameter pair (813, 5), measuring the oscillation probability 
for v e — > and v e — » will generally result in two allowed regions of the parameter 
space. The first one contains the physical input parameter pair and the second, the 
"intrinsic ambiguity", is located elsewhere. Worse than that, new degeneracies have 
later been noticed ED], resulting from our ignorance of the sign of the Am^ tm 
squared mass difference and from the approximate [^ 23 ,7r/2 — # 23 ] symmetry for the 
atmospheric angle. In general, for each physical input pair the measure of P(v e — > z/ M ) 
and P(v e — > Ufj) will result in eight allowed regions of the parameter space (if the sign 
of Am 2 atm and the 623- octant will still be unknown by the time the Neutrino Factory 
will be operational) |3Ti] . 

Different proposals have been suggested to deal with the three ambiguities. The sign 
of Am 2 atm could for example be determined by combining two of the planned superbeam 
facilities, one of them with sufficiently long baseline while the second one with good 
#13 sensitivity [31]. The #23-octant could be determined by combining one of these 
superbeam facilities (e.g. JHF-I) with a reactor-driven detector [32] (although in this 
case systematics could represent a serious issue). Finally, the intrinsic ambiguity can 
be solved by fitting at two different baselines at the same time [28j; or by increasing the 
energy resolution of the detector j33]-[3n]- in [HZ! a detailed study of the combination 
of a superbeam and of a (single-baseline) Neutrino Factory was presented, showing 
that this combination is extremely helpful in solving ambiguities whenever a 1 Mton 
superbeam-driven water Cerenkov detector is added to the Neutrino Factory-driven 40 
Kton magnetized iron detector (MID) that was considered in previous studies [3%] . 

In all of these proposals only the "golden channel" (v e — > v^, or — > u e in the case 
of superbeams), with different detectors and neutrino sources, was considered. How- 
ever, in pP it was noticed that muons proceeding from r decay when r's are produced 
via a v e —>■ v T transition show a different (813,6) correlation from those coming from 

- * *V By using a lead-emulsion detector based on the Emulsion Cloud Chamber 
(ECC) technique, capable of the r-decay vertex recognition, it is therefore possible to 
use the complementarity of the information from u e — > v T and from v e — > to solve 
the intrinsic (9is,S) ambiguity. The lesser statistical significance of the former 1 with 

1 Due to the mass of the ECC detector, to the v T N CC cross-section, to the BR(t — > /i) factor and 
to the need for the vertex identification, which results in a smaller total efficiency. 
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respect to the latter (the "golden" muons), inspired the nickname of "silver channel" 
and consequently of "silver" muons. Silver muons, in combination with golden muons, 
could also be extremely helpful in dealing with the [0 2 3, 7r/2 — #23] ambiguity, since the 
leading term in P(v e —> v T ) is proportional to cos 2 #23, whereas the analogous term in 
P{y e ~^ Vy) is proportional to sin 2 #23. However, the sensitivity of the silver/golden 
channel combination to the ^23-octant strongly depends on the value of #13. A detailed 
study of this combination to solve the ^-octant ambiguity is currently underway 
(sec [39 ). In this paper, as in [T], we restrict ourselves to the (#13, 5) ambiguity, by 
fixing #23 — 45° and by choosing a given sign for Am 2 tm (in the hypothesis that more 
information on the three neutrino spectrum will be available by the time the Neutrino 
Factory will be operational). Clearly, solving the three ambiguities at the same time 
will need the combination of different measurements. 

In jT], a preliminary analysis of the foreseeable backgrounds in the measurement of 
p e — > v T transitions was used to substantiate the proposal with a realistic simulation 
of the detector based on [T3J ED]. The aim of this paper is to present a dedicated 
analysis to measure the "silver" muon signal at an OPERA-like lead-emulsion detector. 
Furthermore, a detailed study of the systematic errors associated to both ECC and MID 
is also presented. In particular, it is shown that, by including a realistic systematic error 
in the analysis of the golden channel, its performance in determining the (6*13, S) value 
is considerably worsened if compared to the ones reported in [38J. A refined statistical 
analysis of "silver" and "golden" muons is also included, in order to optimally deal 
with very low statistics signals but with a very high signal-to-noise ratio. We consider 
in this analysis two possible combinations of detectors and baselines: 

• An ECC at L = 732 km and a MID at L = 3000 km 
(same combination considered in PQ); 

• An ECC at L = 3000 km and a MID at L = 3000 km . 

Our conclusion is that the combination of "golden" and "silver" channel at the 
Neutrino Factory can indeed solve the intrinsic ambiguity problem for # 13 > 1°, thus 
confirming the results of PQ. For 13 = 1° we are on the edge of the experimental 
sensitivity of the silver channel with both configurations. However, given the lower 
expected background, we find that the second configuration (both detectors at 3000 km) 
allows a reduction in the size of the confidence interval for #13 > 1°. 

The paper is organized as follows: in Sect. El we introduce the "golden" and "silver 
channels" at the Neutrino Factory and we show how a combination of these two signals 
could solve the intrinsic ambiguity (#13, 5); in Sect. EH we recall the neutrino fluxes, the 
vN cross-sections and finally report the expected number of events for each channel; 
in Sect. 0] the characteristics of a lead-emulsion OPERA-like detector are depicted; in 
Sect. E] we discuss in full detail the "silver muons" signal and the different sources of 
background at an ECC detector; in Sect.|H]we give an estimate of the required scanning 
load; in Sect. [7| we perform a refined analysis of low statistics signals with a very good 



4 



S/N within the Feldman- Cousins 41J approach; in Sect. [HI we eventually draw our 
conclusions. 

2 The Golden and the Silver Channel 

Following eq. (1) of [28J, we get for the transition probability u e — > (y e — > v^) and 
u e — ► f T (i/ e — > u T ) at second order in perturbation theory in 8i 3 , A Q /A atm , A Q /A and 
A Q L (see also 02]- EH), 

P±(0 13 , 5) = X ± sin 2 (2fl 13 ) + Kt cos(fl 13 ) sin(2fl 13 ) cos (±5 - ^f^j +Z, (2) 

and 

P±(0 13 , 5) = XX sin 2 (2# 13 ) - F± cos(£ 13 ) sin(2# 13 ) cos (±6 - + Z T , (3) 

where ± refers to neutrinos and anti neutrinos, respectively. The coefficients of the 
two equations are: 




with A = \f2Gpn e (expressed in eV 2 /GeV) and B T — \A =p A atm | (with =f referring 
to neutrinos and anti-neutrinos, respectively). Finally, A atm = Am 2 atm /2E V and A = 
Am%/2E V . 

The parameters #i 3 and 5 are the physical parameters that must be reconstructed 
by fitting the experimental data with the theoretical formula for oscillations in matter. 
In what follows, the other parameters have been considered as fixed quantities, sup- 
posed to be known with good precision by the time when the Neutrino Factory will be 
operational. In particular, 

1. in the solar sector we fixed 812 = 33° and Am 2 & = 1.0 x 10~ 4 eV 2 (same as in pQ). 
Although these values do not correspond to the present best-fit value from the 
combined analysis of the solar and KamLand data [TO], they have been chosen 
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to make direct comparison with our previous results. Furthermore, notice that 
the squared mass difference lies in between the two allowed solutions LMA-I with 
Am| = 7 x 10" 5 eV 2 and LMA-II with Am| = 1.4 x 10" 4 eV 2 g]. 

2. in the atmospheric sector, # 2 3 = 45° and Am 2 tm = 2.9 x 10~ 3 eV 2 7 , with Am 2 tm 
positive. Notice that for 6*23 = 45° the 6* 2 3-octant ambiguity [3(J is absent. 

Finally, we also considered a fixed value for the matter parameter, A = 1.1 x 10~ 4 
eV 2 /GeV (this value, obtained by using the average matter density alongside the path 
for the chosen distance computed with the Preliminary Earth Model [12] > is a good 
approximation for the case under consideration of L < 4000 km). For simplicity, we 
have not included errors on these parameters 2 . 

Eqs. (J2J) and © lead to two equiprobability curves in the (6*13, 5) plane for neutrinos 
and anti-neutrinos of a given energy: 

Pi(W)=P e t(W); (7) 



P±(W) = P^(0 13 ,5). (8) 

Notice that X± and Z T differ from the corresponding coefficients for the v e — > 
transition for the cos #23 ^ sin #23 exchange, only (and thus for #23 = 45° we have 
X = X T , Z = Z T ). The Y± term is identical for the two channels, but it appears 
with an opposite sign. This sign difference in the F-term is crucial, as it determines a 
different shape in the (#13, 5) plane for the two sets of equiprobability curves. 

In Fig. Q we superimposed the equiprobability curves for the u e —>■ v T and v e — > 
oscillations at a fixed distance, L = 732 km, with input parameters #13 = 5° and 
5 = 60°, for different values of the energy, E u e [5, 50] GeV. The effect of the different 
sign in front of the F-term in eqs. (J2J) and Q can be seen in the opposite shape in 
the (#i3, 5) plane of the u e — > v T curves with respect to the v e — > ones. Notice that 
both families of curves meet in the "physical" point, #13 = ^13, 5 = 5, and any given 
couple of curves belonging to the same family intersect in a second point that lies in a 
restricted area of the (A#, 5) plane, the specific location of this region depending on the 
input parameters (^13, 5) and on the neutrino energy. Using a single set of experimental 
data (e.g. the golden muons), a x 2 analysis will therefore identify two allowed regions: 
the "physical" one (around the physical value, ^13, 5) and a "clone" solution, spanning 
all the area where a second intersection between any two curves occurs. This is the 
source of the intrinsic ambiguity pointed out in [2H]- When considering at the same 
time experimental data coming from both the golden and the silver muons, however, 
we see that the clone regions for each set of data lie well apart. Thus, a comprehensive 
X 2 analysis of the data will in principle result in the low-x 2 region around the physical 
pair, only. Of course, this statement is only true if the statistical significance of both 
sets of experimental data is sufficiently high. 

2 It has been shown in |2B] that the inclusion of the foreseeable uncertainties on these parameters 
does not modify the results on the #13 and 5 measurements in a significant way. 
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3 Expected rates at the Neutrino Factory 

To get the expected number of events per bin, we must now convolute the neutrino 
fluxes at the Neutrino Factory with the charged-current (CC) cross-section and with the 
transition probability in eqs. (J2J) and Q. At low energies the neutrino scattering cross- 
section is dominated by quasi-elastic scattering and resonance production (see [HH H?j 
and references therein). However, if E v is greater than ~ 10 GeV, which is the case for 
high-energy muons in the storage ring (30-50 GeV), the total cross-section is dominated 
by deep inelastic scattering (DIS) and can be approximately described with the cross- 
section on an isoscalar target [HJ : 

<r(v + N^r+X) 
a(P + N + X) 

Eq. (JOJ) is not a good description of the DIS cross-section for v T interactions in the 
considered range of energies, since the lepton mass cannot be neglected. Taking into 
account these effects, the total cross-section for the u x (v x ) — N interactions appears 
as in Fig. |21 where we show its behaviour as a function of the neutrino energy. The 
theoretical number of expected charged-current events (with production of final leptons 
£p) in the detector, for a fixed parent muon energy E^, is: 

N tp (y a -> v p ) = N Kton xN A x 10 9 ^ dE v a Vg {E v ) P a p(E v - 6 13 , 5) (10) 
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Figure 2: Behaviour of the (anti)neutrino-nucleon cross section as a function of the 
neutrino energy for z/ e (z/ e ) (long dashed lines), v^ip^) (continuum lines) and z/ T (z/ r ) 
(dashed lines) 

where Nxton is the detector mass in Kton and Na x 10 9 is the number of nucleons per 
Kton. 

In Tab. Q we report the neutrino and anti-neutrino charged- currents interaction 
rates assuming no oscillation for two different baselines (namely L = 732 and 3000 km) 
for an ideal (i.e. with perfect efficiency) detector with a mass of 1 Kton. We consider 
three different muon polarization values V^=0, ±0.3 (the "natural" polarization) and 
±1. The neutrino beam results from the decay of 2 x 10 20 yU +, s and /i~'s per year in a 
straight section of an = 50 GeV muon accumulator, with five operational years for 
each polarity. Fluxes have been integrated in the forward direction with an angular 
divergence (taken to be constant) Sip ~ 0.1 mr. Beam divergence effects and QED one- 
loop radiative corrections to the neutrino fluxes have been taken into account in [BP] . 
The overall correction to the neutrino flux has been shown to be very small, O(0.1%). 

The situation in the hypothesis of neutrino oscillations at the same baselines is 
given in Tabs. 121 and El The considered PMNS matrix parameters are: 12 = 33°, 
Am 2 12 = 10~ 4 eV 2 ; 6 23 = 45°, Am 2 3 = 2.9 x 10~ 3 eV 2 ; 9 X3 = 5° and 5 = 90°. 

4 Detector issues 

In this Section we discuss the basic performance of an OPERA-like detector at the 
Neutrino Factory in reconstructing neutrino interactions. In doing so, we profit of the 
experience and of the test results gathered within the OPERA Collaboration HU] • 
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Table 1: Neutrino and anti-neutrino charged- currents interaction rates for L = 732 km 
and 3000 km per 1 Kton and per 5 operational years when 2 x 10 20 muons decay 
in the straight section of the storage ring. These fluxes have been averaged over an 
angular divergence of 0.1 mr. The results are presented for three muon polarizations 
V^= 0, ±0.3 and ± 1. 
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Table 2: Oscillated charged- currents event rates for fi~ (upper table) and fi + (lower 
table) beam assuming neutrino oscillations with 6>i 3 = 5° and 5 = 90° in a 1 Kton 
detector, for a L = 732 km baseline for different polarizations of the parent muon. We 
have considered 1 x 10 21 muons decays (2 x 10 20 useful muons /year x 5 operational 
years). 
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Table 3: The same as Tab. but for L = 3000 km. 



4.1 The detector layout 

The main features of the detector we are going to discuss in this paper are very similar 
to the ones of the OPERA experiment fHHHI], which is meant for a long baseline search 
for —>■ u T oscillations at the CNGS beam. The experiment uses nuclear emulsions 
as high resolution tracking devices for the direct detection of the r produced in the v T 
CC interactions with the target. 

OPERA is designed starting from the Emulsion Cloud Chamber concept (see refer- 
ences quoted in ^3|) which combines the high precision tracking capabilities of nuclear 
emulsions and the large mass achievable by employing metal plates as a target. The 
basic element of the ECC is the cell which is made of a 1 mm thick lead plate followed 
by a thin emulsion film. The film is made up of a pair of emulsion layers 50 /im thick on 
either side of a 200 fim plastic base. Charged particles give two track segments in each 
emulsion film. The number of grain hits in about 50 /xm (15-20) ensures redundancy in 
the measurement of particle trajectories. By piling-up a series of cells in a sandwich-like 
structure bricks can be built, which constitute the detector element for the assembly 
of massive planar structures (walls). A wall and its related electronic tracker planes 
constitute a module. A supermodule is made of a target section, which is a sequence of 
modules, and of a downstream muon spectrometer. The detector consists of a sequence 
of supermodules (see Fig. EJ). 

The signal of the occurrence of a v T CC interaction in the detector target is identified 
by the detection of the r lepton in the final state through the direct observation of its 
decay topology. A r may decay either into the lead plate where the interaction occurred 
(short decay) or further downstream (long decay). For long decays, the r is detected 
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Figure 3: OPERA at CNGS. 



by measuring the angle between the charged decay daughter and the parent direction. 
The directions of the tracks before and after the kink are reconstructed in space by 
means of the pair of emulsion films sandwiching the lead plate where the interaction 
occurred. A fraction of the short decays is detectable by measuring a significant impact 
parameter (IP) of the daughter track with respect to the tracks originating from the 
primary vertex. 

The detection of the r decay and the background reduction benefit from the dense 
brick structure given by the ECC, which allows the electron identification through its 
showering, the pion and muon separation by the dE/dx measurement method, and the 
determination of the momentum of each charged particle employing techniques based 
on the Multiple Coulomb Scattering. All these methods are discussed in the following. 

Electronic detectors placed downstream of each emulsion brick wall are used to 
select the brick where the interaction took place (to be removed for the analysis) and 
to guide the emulsion scanning. The target electronic detectors are also used to sample 
the energy of hadronic showers and to reconstruct and identify penetrating tracks. 

In the following we assume a supermodule structure with 31 brick walls followed 
by a muon spectrometer. For the present analysis, we consider the same setup as for 
OPERA at CNGS. The electronic detectors are 1 cm thick, 2.6 cm wide scintillator 
bars. The active detectors at the spectrometers are Resistive Plate Chambers (RPC) 
with 3 cm wide readout strips. The bending of the particle trajectory before and after 
the magnet is measured by high precision drift tubes (DT) ■ With this configuration 
the total target mass of a supermodule is 0.9 Kton. 
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For details on the event reconstruction both with the nuclear emulsions and the 
electronic detector, and on the OPERA sensitivity to — > v T and — > v e oscillations, 
we refer to [3 El HO] ■ 



4.2 Kinematical analysis with ECC 

4.2.1 Momentum measurement through the multiple scattering 

When a particle of momentum p and velocity (3 traverses a material of thickness X, 
measured in units of radiation length X , the distribution of the scattering angle 6q in 
a plane is approximately gaussian with a RMS given by 

One can determine 8q by measuring the distribution of the angular difference in two 
consecutive emulsion films. This method is called angular method. A resolution 89 of 
about 2 to 3 mrad can be routinely obtained as shown in Fig. HJ The angular resolution 
is affected by both systematical sources (i.e. planarity of the scanning surface) and 
statistical effects which can be significantly reduced by means of precise and repeated 
measurements respectively providing a resolution of about 1 mrad. 

The RMS of the measured scattering angles, 6 2 M) is the quadratic sum of the scat- 
tering signal #o and of the measurement error SO: 9 2 M = 8q + 58 2 . The error on 9q when 
iV independent measurements are performed is 



59 n 



(12) 



Using Nfu m emulsion films, N = (Nfu m — 1) x 2 points can be used, since each 
of the two projections provides a set of independent measurements. Therefore, one 
can easily obtain a Sp/p ~ 10% with a brick of about 50 emulsion films. In Fig. we 
report the momentum resolution corresponding to an angular resolution of 3 mrad for 
2, 3 and 4 GeV negative pions impinging onto a 3Xo ECC (17 emulsion films) jSI]. A 
resolution of 28% (2 GeV) to 36% (4 GeV) has been achieved. 

4.2.2 Electron identification 

A typical example of an electron identified in an ECC brick is shown in Fig. El The 
identification is done employing the multiple scattering technique discussed above ap- 
plied to the electron track before showering and counting the number of tracks when 
the shower develops (calorimetric measurement). In particular the peculiar energy loss 
rate of electrons is used to discriminate between electrons and hadrons: the energy 
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Figure 4: The measured angular resolution routinely achieved with an ECC shown for 
both angular projections. 



is practically unchanged for hadrons but it is sizebly decreased by bremsstrahlung for 
electrons while traversing the brick. A x 2 estimator to test the two possible hypothe- 
ses (hadrons or electrons) is built and it is used as a separator. The \ 2 minimization 
provides also the particle energy estimation. The other way to identify electrons is to 
count segments associated to its shower. The lateral spread of a shower is of the order 
of 1 cm which makes the shower itself well confined within a single brick. An energy 
resolution of about 20% can be obtained with this method, limited by background 
tracks falling into the analysis region. 

Preliminary studies show that combining the two methods it is possible to achieve 
97% identification efficiency for electrons with energy above 1 GeV |13| . 
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Figure 5: T/ie 1/p resolution achieved with a 3X ECC and 3 mrad angular resolution 
exposed to 2, 3 and 4 GeV pions. 



4.2.3 Low momentum muon identification 

Low momentum muons stopping in the ECC can be identified measuring the energy 
loss rate near the end-point of the range In this regime, the energy loss rate is 
given by I = dE/dx = k//3 2 . A measurement of the energy loss rate I and of the 
distance x traveled by the particle before stopping will therefore give the particle mass 
M = (2/j3 2 )E = (2/k)I 2 x, with an error that gets contributions from 5x and 51. The 
latter is the dominant one when a particle stops in the brick. 

The grain measurement amounts to about (30 ± v / 30)/100 /im; hence, by mea- 
suring only one emulsion sheet one can get about 18% resolution on the ionization 
measurement. Since muons and pions are expected to cross about 10 emulsion films 
after entering into the non-relativistic regime, a mass resolution 5M/M ~ 251 /I = 
2 x (0.18/\/l0) ~ 0.12 is achievable which means 5M = 16 MeV for a pion. This mass 
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Figure 6: Electromagnetic shower observed in a brick exposed to an electron beam. 



resolution is used to separate muons from pions. 
4.3 Muon reconstruction 

The magnetic spectrometers located in the back of the ECC target are aimed to recon- 
struct the charge and the momentum of penetrating particles. The fraction of "silver" 
candidates where it is possible to measure the muon sign and the corresponding mis- 
assignment probability is a key parameter to assess the performance of OPERA as 
a detector for the Neutrino Factories. In fact the primary muon identification, even 
without charge reconstruction, plays an important role to veto background events. In 
particular, charm-production from CC interactions can be significantly suppressed 
if a prompt muon from the primary vertex is identified. Both at CNGS and in the 
present analysis muon candidates are selected combining the information of the scin- 
tillator trackers and the spectrometer RPC's through a pattern recognition algorithm. 
For a detailed description of this procedure and the build-up of the candidate tracks we 
refer to ^3] (Section 5.4) and j3U] (Section 5.3). The muon identification algorithm is 
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described in details in [T3] (Section. 7.2) and [10] (Section 5.4) and it is briefly recalled 
here. 

4.3.1 Muon identification 

Among the reconstructed tracks, the longest one is chosen as a possible muon candidate. 
If the track is very penetrating and exits from the back of the spectrometer, it is 
validated as muon and no further cuts are applied. Otherwise cuts are imposed on 
the number of target walls and spectrometer iron plates belonging to the candidate 
and on the isolation of the track. Finally, a candidate is validated if the momentum 
measured by range is consistent (within proper tolerances) to the one measured in 
the ECC through Multiple Coulomb Scattering. All the ECC tracks matching the 
muon candidate within 200 mrad are considered. If more than one ECC track exhibits 
range-momentum correlation, the one with the best angular match is selected. The 
performance of the algorithm can be expressed in term of muon identification efficiency 
and matching probability. The former represents the number of true muons passing 
the /i-id cuts and matched to a track in the emulsions. The latter quantifies the 
probability of correct matching to ECC. Fig.[7|shows the muon identification efficiency 
and the matching probability versus the energy of the muon in v T — ► t~X — > fi~v^v T X 
interactions. It is worth noting that the muon identification efficiency for short muons 
stopping before the spectrometer can be improved significantly measuring the dE/dx 
of the particles in the emulsions around the range-out area (see Section f4. 2. This is 
particularly interesting for charm vetoing and in the OPERA experiment running at 
the CNGS allows a reduction of the charm background by about a factor 2. However, 
these results are only based on Monte Carlo simulations and tests are in progress in 
order to measure experimentally the ir/fj, separation through the dE/dx technique. 
Although preliminary results are in a good agreement with Monte Carlo estimates, 
conservatively for the present analysis this technique is not employed. However, as an 
example, we will give results by using this technique only for a given (#i 3 , 8) pair, see 
Section 1731 . 

4.3.2 Charge measurement 

OPERA is able to measure the charge of the muons crossing the magnetic spectrometers 
or stopping in them. For through going muons an iterative \ 2 minimization procedure is 
used to fit the measured muon trajectory combining the hits from target scintillators, 
spectrometer RPC's and drift tubes. Hence, both the particle momentum and the 
charge can be reconstructed. For stopping particles the best estimate of the momentum 
is provided by the range, while the charge reconstruction still relies on the curvature in 
the magnetic field. The fraction of events where the charge can be measured (through 
going muons and stopping muons with enough spectrometer hits) and the probability 
to assign correctly the charge are shown in Fig. 0as a function of the muon energy. As 
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already discussed in [^] the charge mis-assignment probability in the energy range of 
interest is approximately constant and it is at the level of 0.1-0.3%. 
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Figure 7: Percentage of v T — > t~X — ► /i~z/ M z/ r X events with muon identified (dots), 
well matched to the ECC track (square), reaching the spectrometer (triangle) and with 
charge correctly assigned (crosses) versus the energy of the muon. To ease the reading, 
the width of the bins and the statistical errors are shown just for the first mark. 



5 Analysis of silver muons 

Silver wrong-sign muons are produced in z/ T CC interactions, coming from v e oscillations 
with subsequent muonic decay of the r. In the following we discuss in detail the signal 
efficiency and the expected background under the hypothesis that /z + 's circulate into the 
muon storage ring (namely, a 50% v e plus 50% beam hits the detector). Therefore, 
the signal we are looking for is v e — >• v T — > t~ — > . 
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5.1 Expected number of signal events 

The search for muonic tau decays is performed by using an approach very similar to the 
one described in Refs. |13| I4U|. Two classes of events are considered: short (the r decays 
into the lead plate where the interaction occurred) and long decays (the decay occurrs 
outside the vertex plate, therefore the kink angle in space can be reconstructed). 



a) b) c) 




Figure 8: Sketch of signal a) and background b) and c) topologies for r — > \i decays. 

The short decay search exploits the impact parameter technique, searching for a 
muon with a large IP with respect to a reconstructed vertex of at least two tracks (see 
Fig.IHJ). The main background is a subset of charm-production. Charm can be produced 
in CC and contributes to the signal either when the charge of the primary muon 
is wrongly measured in case of neutral charmed-hadron production (Fig. IBl-b) or when 
the primary muon is not identified in the case of charged charmed-hadron production 
(Fig. IBl-c). Charm can also be produced in (from v e — > v^) CC interactions. If the 
origin of the muon track is unknown or wrongly reconstructed, these events can mimick 
a r — > n decay (as sketched in Fig. [SJ). 

This background is significantly reduced by imposing a lower cut at 2 GeV on the 
invariant mass of the hadronic system, as computed by the momenta measured in the 
ECC. Even with a modest momentum resolution of 50% the charm background can be 
reduced by more than a factor of 1000 retaining 15% of the r — > \i short decay signal. 
The results reported in show that the total efficiency for muonic short decays is 
about 0.7%, including the muonic branching ratio of the r. Although this analysis is 
very preliminary and could be improved, in the following we assume such an efficiency. 
The expected background was evaluated for the CNGS conditions to be 2 x 10~ 6 x ^ M CC 
DIS events, which is very small. Moreover, having in mind that at a Neutrino Factory 
the bulk of the events comes from P M interactions, while at the CNGS from u^, we can 
safely assume that the short decay channel is background free. 

For long decays, r candidates are searched for by measuring the kink angle in space 
and taking into account the worse angular resolution for decays in the base. The in- 
efficiency in the kink detection results from the need of rejecting small-angle and very 
large-angle kinks. The upper cut (> 500 mrad) is motivated by considerations related to 
the scanning efficiency of the automatic microscopes. The lower cut (20 mrad in space) 
has been assigned taking into account that experimental data for hadron rescattering 
are available only down to such an angle. The excellent intrinsic angular resolution of 
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the emulsion films would allow a substantially looser cut. Nevertheless, awaiting back- 
ground measurements extended to lower angles, we maintain the conservative 20 mrad 
cut. As in the analysis described in Ref. ^3] we only consider decays occurring in two 
lead plates downstream from the vertex plate. Unlike the CNGS analysis described 
in the OPERA Proposal, we do not apply neither a cut on the muon momentum nor 
a cut on the visible energy. On the other hand the charge measurement of the muon 
track is required. The only kinematical cut we apply is the transverse momentum at 
the decay vertex to be larger than 250 MeV. This cut allows the rejection of muonic 
decays of pions and kaons. 

Summing up the short and long decay contributions, the overall efficiency (including 
the muonic tau decay branching ratio) for an average neutrino energy of 35 GeV is 
about 5%. 

The expected number of signal events detected during a 5 year data taking with a 
1 Kton detector for different values of (6*13, 5) is given in Tab. |U 



(013,*) 


T -► yU 
L = 732 km 


T -> H 

L = 3000 km 


(l°,-90°) 


0.64 


0.72 


(5°, -90°) 


8.73 


8.90 


(i°,o°) 


0.0092 


0.070 


(5°,0°) 


5.57 


5.66 


(1°,90°) 


0.52 


0.30 


(5°, 90°) 


8.11 


6.81 



Table 4: Expected number of muonic tau decays detected into a 1 Kton OPERA-like 
detector for various (#13, 8) values and baselines (i.e., the t^ 1 — > fi^ branching ratio and 
the overall estimated efficiency are taken into account). We have considered 1 x 10 21 
muons decays with "natural" polarization (2 x 10 20 useful muons/year x 5 operational 
years) for each polarity. 



5.1.1 The golden muons signal at an ECC detector 

Although not competitive with large magnetized calorimeters, an OPERA-like detector 
can also be used to study u e — > oscillations. From an experimental point of view 
this search does not imply additional efforts with respect to the ones needed to study 
the silver channel. Indeed, as soon as a wrong sign muon is reconstructed in the 
detector the brick where the interaction occurred is removed and the event carefully 
analyzed. If a decay candidate is found the event falls in the silver channel sample and 
the appropriate kinematical analysis applied. Conversely, it is classified as a golden 
event. 
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In an OPERA-like detector the background from charmed hadron and r (produced 
by v T coming from v e — > v T oscillations) decays is highly suppressed because of the 
detector capability in detecting decay topologies. Since the presence of a kink on the 
muon track is not required, the background from h~ decays or h~ punch-through in 
this channel is higher. 

Given the fact that, also in the most optimistic case, the mass of the OPERA-like 
detector is small if compared to a MID detector and it is not fully magnetized, the 
gain in sensitivity one could have in studying this channel is modest. Therefore, in the 
following we focus only on the v e — > v T sensitivity and we use it in combination with a 
MID detector exploiting the golden channel. 

5.2 Expected background in the long decay sample 
5.2.1 Neutrino induced charm-production 

Charm-production from neutrinos can be induced by both unoscillated v e and com- 
ing from oscillated v e through the reaction 

uiN r + C + + X 

If the primary lepton is not identified and the charge of the muon from the semi- 
leptonic decay of the charmed hadron is wrongly measured, then the event is classified 
as a silver candidate. 

The charm-production rate has been estimated by using the results of a combined 
analysis of all available data on neutrino- induced charm- production [53] • By using the 
z/^, v e predicted spectrum, we expect a charm rate, normalized to CC interactions, of 

R c {y^ = (5.12 ± 0.30)%; R c (u e ) = (4.74 ± 0.23)% . 

Given the capability in detecting decay topologies, an ECC detector is only sensitive 
to semi-muonic decays of the charged charmed mesons, unlike the electronic detectors. 
As it can be derived from Ref. [51], the fraction of charged among produced charmed 
particles is, at the Neutrino Factory energies, of about 45% (fc+)- Moreover, given 
the fact that the search is limited to muonic tau decays, we have to consider only the 
semi-muonic branching ratio of charmed hadrons, which amounts to about 10% |55j . 

Finally, we can write the expected number of events from this background as 

N c = Nf c x R c (l) x f c+ x (1 - euo) x BR(C + -> /+) x e^ +ID x (1 - e charge ) x e det 

where Nf c is the total number of charged-current events induced by vi, eun is 
the efficiency to identify the primary lepton, e^+jc is the probability to identify the 



20 



muon produced at the decay vertex and to measure the charge for muons reaching 
the magnet spectrometer, e c h arge is the probability to correctly identify the charge of 
the daughter muon and e^et is the probability to detect the charm decay. The latter 
efficiency includes the trigger efficiency, a fiducial volume cut (emulsions cannot be 
scanned up to the edge and there are small cracks between bricks, see Ref. [T3| for 
more details), the probability to correctly identify with electronic detectors the brick 
where the interaction occurred, the efficiency of the emulsion tracking algorithm in 
reconstructing the interaction vertex in the ECC. 

In our case only prompt unoscillated u e contribute to the background, although this 
background is highly suppressed by the fact that the charge of the muon produced at 
the charm decay vertex is positive. Therefore, it contributes to the background only if 
the charge of the muon is wrongly measured. By assuming an electron identification 
efficiency of 97% and the muon identification and charge determination efficiencies 
coming from the algorithm described in Section \4.'S\ we expect a background smaller 
than 10" 8 x Nf c . 

5.2.2 Anti-neutrino induced charm-production 

Charm-production from anti-neutrinos can be induced from both unoscillated and 
v e from — > v e oscillations through the reaction 

UiN -> l + + C~ + X . 

It is worth noticing that this background does not profit of the charge measure- 
ment of the muon from charmed-hadron decay, being the charge of the decay product 
negative. 

The charm-production rate in anti-neutrino induced charm-production has been es- 
timated by using the approach discussed in Appendix A. The expected rate, normalized 
to CC interactions, are 

Bjy^ = (4.84 ± 1.94)%; R c (v e ) = (4.36 ± 1.48)% . 
Finally, we can write the expected number of events from this background as 

- rr - - - ' ~ ' 

N c = iV7° x R c (l) x f c + x BR(C -»• I ) x (1 - s lID ) x e charm x e^- ID xe det 

where gives the product of the probability that the primary lepton is not identi- 
fied, the decay is classified as long and the secondary muon is reconstructed with the 
right charge. This probability has been evaluated to be 1.7 x 1CT 3 and it is mildly 
dependent on the neutrino energy. 
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In our case only prompt unoscillated contribute to the background, which amounts 

to 

N c = 3.7 x 10~ 6 x Nf c . 

5.2.3 Events from r + — > /i + decays 

A potential source of background comes from r + , produced by v T from the leading 
Vfj, — > v T oscillations and undergoing to a charged- current interaction, decaying into a 
positive muon with the charge wrongly measured. 

From Tables 121 and El we see that for 6 13 = 5° and 5 = 90° the number of r + is 
one order of magnitude larger than t~ . Therefore, under the assumption that the r- 
decay detection efficiency is independent of the lepton charge and that the probability 
to misidentify the charge of the muon is of the order of 10~ 3 , we expect that for 
this particular value of (#13, 5) the background from this channel is a factor hundred 
smaller than the observed number of signal events. Nevertheless, the dependence on 
the oscillation parameters is properly taken into account. 

5.2.4 Muons wrongly matched to a hadron track 

Occasionally, a fi + identified as /i~ in the electronic detector is wrongly matched to a 
hadron track at the exit of the vertex brick, as it is reconstructed by tracking with the 
emulsion ECC brick. Because the kinematical analysis in the muonic channel is very 
loose, without special care in the track matching such a mismatch would result in a 
non negligible number of background events due to hadron re-interactions. 

In order to reduce this background, if there is an emulsion track with an angular 
difference smaller than 50 mrad with respect to the one matched to the "electronic" 
muon, the matching is flagged as ambiguous and, in the presence of a kink, the event is 
not taken as a r candidate. This implies a few percent relative loss of efficiency already 
accounted for in Section I5~T1 Further improvements in the background reduction come 
from the requirement to have the charge of the muon correctly identified, the latter 
being not present in the OPERA analysis discussed in Ref. |4T)] . 

The mismatch probability has been computed to be 6 x 10 -4 . On the other hand 
the probability that a hadron undergoes to a re-interaction mimicking a kink topology 
which survives the kinematical cuts has been computed with the FLUKA package jSHl 
to be 8 x 10~ 4 . Therefore, the expected background from this source is 

N mism = 7 x 10- 9 x . 
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5.2.5 Decay in flight of h and punch-through h 

Neutral-current (NC) events with punch-through hadrons or particles decaying in flight 
will dominate the scanning load of an OPERA-like detector at the Neutrino Factory 
(see Section EJ). Moreover, if the hadron undergoes to an elastic scattering in lead 
mimicking a kink topology, it will also contribute to the silver channel background. The 
probability of a primary hadron to be reconstructed as a muon, reach the spectrometer 
and have a reconstructed charge consistent with a silver candidate has been computed 
by Monte Carlo simulation. It corresponds to an efficiency of 0.6%. Folding the events 
with the probability of undergoing large px scattering in the first two lead sheets we 
obtain an efficiency of 4.1 x 10~ 6 . About 83% of these events are decay in flight. The 
other 17% fraction is made up of punch-through pions surviving the cut on the muon 
length. Also u e CC events contribute to this background if the primary electron is 
unidentified. The background is 

Ndecay = 4.1 X 1(T 6 X E det X 

where Nf c (Nf c ) is the expected NC (CC) rate for neutrinos of type I; e e is the 
electron efficiency and Edet is defined as in Section 15.2.11 

5.2.6 Large-angle muon scattering 

A muon can undergo to a large-angle scattering mimicking a muonic decay of a short- 
lived particle. Potential sources of this background are produced in CC interactions 
of from v e — ► oscillations and /i + from z^CC interactions, with a wrongly mea- 
sured charge. 

Extensive studies of this background have been carried out by the OPERA Collab- 
oration. A Monte Carlo simulation including the lead form-factors gave a rate of muon 
scattering off 2 mm lead mimicking a r decay of 0.2 x lO^/N^ [13]. Recently, the 
analysis of a dedicated measurement of large-angle scattering of 9 GeV muons in lead 
plates became available. The corresponding measured rate is (0.6lo g) x 10~ 5 /N^ c |4H] . 
In the OPERA proposal, a conservative assumption of 1 x 10~ 5 /N^' C was made. Even 
in this case, the background from this source is smaller than 0(1O~ 8 ) x N^ c 

5.2.7 Associated charm-production 

Another possible source of background is given by the associated charm-production, 
in which two charmed hadrons are produced and one of the two escapes the detection. 
The cross-section of this process has been measured only by the E531 experiment and 
it turns out to be more than one order of magnitude smaller than that for single 
charm-production. In the following, we neglect this background. 
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6 Estimate of the scanning load 



We try here to give an estimate of the overall scanning load needed to handle the silver 
and golden muon signals at the ECC detector and how it scales with the detector mass. 

There are several sources of events contributing to the scanning load. Considering 
a fi + in the storage ring, we can classify them in four categories: signal events from 
Ve — * Vft, v t oscillations; background events from — > z/ M , v T — > with misidentified 
muon charge (as discussed in Sections 15 . 2 . 31 and 15 . 2 . 6|) ; background events from both 
v e and v ^ with charged charmed meson production and unidentified primary lepton (as 
in Sections 15.2.11 and I5.2.2|) ; background events from v e and NC interactions with 
punch-through mesons misidentified as muons (Section I5.2.5j) . Whereas the signal 
events depend on the mixing matrix parameters, the different sources of background 
are mainly parameter-independent and can be estimated by the knowledge of the v e 
and Vfj, fluxes at a given distance from the source. 

Consider a fi + in the storage ring and the mixing matrix parameters of Tab. 12 
namely #13 = 5° and 5 = 90°, a baseline L = 732 km and the "natural" polarization 
Vfj,+ = 0. Any time a (a particle with the characteristics of a muon and with 

charge identified as negative) is reconstructed by the electronic detector, we will have 
to scan the emulsions of the brick where vN interactions occurred. The first category of 
events are signal events: for a 1 Kton detector, we expect ~ 310 events from v e — > 
oscillations and ~ 30 events from v e — > v T oscillations with subsequent muonic tau 
decay, for the considered mixing matrix parameters. In the second category we find 
~ 2200 events from — » with misidentification of the muon charge (however, after 
the emulsion scanning only large angle muon scattering events will contribute to the 
background, Section I5.2.6|) and ~ 2 events from —>■ v T . In the third category we 
find ~ 80 events from v^N — > fi + C~ with non-observation of the primary right-sign 
muon and correctly identified secondary wrong-sign muon from the C~ decay. All 
other kind of events in this category [y e — ► l~C + and — > e + ,r + C~) give negligible 
contributions to the scanning load due to the oscillation probability or to the 3 x 10~ 3 
factor for charge misidentification. These events will be scanned and something like 95 
% of them will be rejected when the primary lepton is identified during the emulsion 
analysis. Finally, events in the fourth category give the most relevant contribution to 
the scanning load: these are NC interactions of both v e and where a pion or kaon 
is identified as a muon (punch-through hadrons or their muonic decay). Considering 
the neutrino fluxes in Tab. Qfor the polarization = and L = 732 km, a factor 
0.3 for the NC/CC cross-section ratio and the probability of 0.6 x 10~ 2 for this process 
to occur (see Section I5.2.5|) we have about ~ 4800 events to be scanned. Given the 
absence of a physical kink these events will be rejected after emulsion analysis except 
for hadron re-interactions. 

In summary, for a 1 Kton mass detector located at L = 732 km from the neutrino 
source, we estimate a total scanning load in 5 years of data taking of less than 1 x 10 4 . 
Notice that the two dominant sources of scanning load, i.e. non-oscillated CC 



24 



interactions with misidentified muon charge and u e , NC interactions with punch- 
through mesons identified as muons, do not depend (to a great extent) on the mixing 
matrix parameters and therefore do not scale with 6 13 . For a 5 Kton detector we 
expect less than 4 x 10 4 bricks to be scanned in 5 years, a scanning load that seems 
technologically affordable. 

Finally, we remind that if the detector were to be located at L = 3000 km from the 
neutrino source, the scanning load will be reduced by about a factor 15. 

7 Sensitivity to (#13, 5) 

The overall sensitivity achievable in the plane (#13, S) is evaluated by combining golden 
muon events measured by a 40 Kton iron detector run with both beam polarities and 
silver muon events observed with a 5 Kton OPERA-like detector. The golden muon 
sample at the OPERA-like detector is not considered in the present analysis since it 
has a worse signal/background ratio compared with the iron detector. Eventually, only 
the /i + polarity is taken into account for the OPERA-like detector. In the following, 
the magnetized iron detector is located at a fixed baseline L = 3000 km, while both 
L = 732 and L = 3000 km are considered as possible baselines for the ECC detector. 

7.1 Signal and background uncertainties 

The total number of golden muons and background events expected for the iron detector 
is reported in (3K|, whereas Section El describes our analysis of the silver muons and 
backgrounds in the OPERA-like detector. In the present section we summarize the 
corresponding systematic uncertainties for both detectors. Although our estimates 
are based upon the present experimental knowledge, whenever a definite experimental 
program allowing an improvement is foreseen we will make use of the expected results. 

The main sources of systematic uncertainties for the "silver" signal events in the 
OPERA-like detector are the knowledge of the emulsion scanning efficiency and the 
cross-section ratio a T /a^. In the following we assume an overall systematic uncertainty 
of 15% which is consistent with the one used in Ref. |13j . 

For an OPERA-like detector three main background contributions (see Section EJ) 
are present: the muonic decay of r + events from — > v T oscillations, the anti- neutrino 
induced charm-production and the decay in flight of h~ and punch-through h~. The 
present knowledge on anti-neutrino induced charm-production, as discussed in Ap- 
pendix A, is limited to about 40% because of lack of data. However, at the Neutrino 
Factory a specific short-baseline program is foreseen jJT], collecting a sample of (9(10 6 ) 
events induced by neutrinos and anti-neutrinos with a charmed hadron in the final 
state. Such a statistics will improve considerably the knowledge of the related cross- 
sections. In the following we assume a conservative 10% systematic uncertainty on 
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the charm background. The uncertainty on the background induced by r + decays is 
the same as the one on the "silver" signal. The last background contribution is dom- 
inated by the poor knowledge of the hadronization and of the hadronic re-interaction 
processes. We conservatively assume a 50% uncertainty on this background source. 

A preliminary study of the performance of a magnetized iron detector for the u e — > 
Ufj, search at a Neutrino Factory was presented in j3Sj. The main background sources 
were pion and kaon decays and muonic decay of charmed hadrons. A 20% resolution 
was assumed on the reconstructed neutrino energy, while no systematic uncertainty 
was assigned to backgrounds and signal events in the calculations. In order to obtain 
realistic estimates of the experimental performance we try to include such uncertainties 
in the present work. The systematic uncertainty on the expected number of signal 
events is mainly related to the detection efficiency and to the calibration, since the flux 
and the cross-section are known to a high accuracy. In the following we assume a 10% 
uncertainty on the signal, taking into account the coarse granularity of the detector 
and the absence of an in situ calibration beam line (used for instance to calibrate 
and monitor the NuTeV calorimeter [SB])- The dependence of the sensitivity on this 
uncertainty will be discussed in the following sections. As far as the background is 
concerned, we assume a 50% systematic uncertainty, consistently with the treatment 
of the OPERA-like detector, although the overall sensitivity only mildly depends on 
this parameter. 

7.2 Statistical treatment 

Since signal and backgrounds have a different energy dependence (see Fig. HJ), the 
number of expected events in both detectors is divided into several energy bins. For 
the iron detector we use five bins of constant width of 10 GeV |3*H1 12 7| l2*Sj. An energy 
resolution of AE/E = 20% is used to smear the expected number of events in each 
bin, resulting in a slightly lower signal/background ratio. The number of signal and 
background events for 5 years data taking and for both polarities in all bins is shown 
in Tab. [5] 

The choice of the binning for the OPERA-like detector is performed according to 
the overall sensitivity to (6i 3 ,5). In particular, since the number of expected events 
at small 9i3 varies quadratically with #13, the binning is optimized for #13 < 2° which 
is the critical region where few signal events are expected and assuming for OPERA 
L = 732 km. Several configurations for both the total number of bins and their 
boundaries are considered. For each binning definition, we evaluate the sensitivity as 
the average 90% CL upper limit obtained by a large ensemble of experiments in the 
absence of a signal. The calculation is performed for different values of 5 in the range 
I 8 |< 60°. The final configuration includes four bins and is the one providing the best 
sensitivity as a function of 5. Fig. fTTH-a shows the sensitivity obtained by changing 
the bin boundaries and Fig. ITUl-b the results obtained with a different number of bins. 
Tab. El summarizes the expected number of signal and background events in each of the 
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Figure 9: Signal (black circles) and backgrounds for a 1 Kton detector as a function 
of the neutrino energy at two different baselines: L = 730 km (top panel), L = 3000 
km (bottom panel). Background from anti-neutrino induced charm-production (white 
circles), r + — > fi + decays (white triangles) and decay in flight of hadrons in and u e 
(white squares) neutral- currents interactions are also shown. The signal corresponds 
to 0i, = 5° and 5 = 90°. 
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Table 5: Number of background and golden muon signal events expected in the iron 
detector at L = 3000 km including our estimate of the systematic effects. The binning 
definition refers to the reconstructed energy . 

reconstructed energy bins. The energy smearing has a larger impact for the OPERA- 
like detector due to the strong energy- dependence of the background (dominated by 
charm-production), thus increasing the expected background at low energy. From 
Tab. El it can be seen that the OPERA-like detector has no sensitivity when both 
#13 < 1° and I 5 |< 15°. This region partially corresponds to the central part of the 
curves in Fig. [THl 
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Table 6: Number of background and "silver" signal events expected in the OPERA-like 
detector at 732 km and 3000 km. The binning definition refers to the reconstructed 
energy. 

In order to evaluate the confidence regions which can be realistically deduced from 
the experimental apparatus, we simulate several sets of data corresponding to different 
(#13; points. For each theoretical point and each energy bin, we throw the numbers of 
(observed) signal and background events from Poisson distributions with mean values 
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Figure 10: Sensitivity (arbitrary units) of the OPERA-like detector as a function of 5 
for #13 < 2° and different bin configurations. The solid curve is the final choice with 
four bins. The dashed lines are obtained by changing the bin boundaries (top panel) 
and by increasing the total number of bins (bottom panel). Signal and background 
uncertainties are taken into account in the calculation. 

corresponding to the expected ones. Furthermore, the mean Poisson values are smeared 
assuming a gaussian width corresponding to the quoted uncertainties (see Section ing . 
Notice that, in principle, the procedure should be repeted many times to compute the 
average confidence intervals from a large ensemble of identical experiments. However, 
given the prohibitive computational time, we only perform two trials for each simulated 
point. 

The individual measurements from each of the 14 energy bins (4 bins in the OPERA- 
like detector for fi + and 5 bins in the iron detector for each beam polarity) are then 
combined using the frequentist approach of Ref. |H] to set 68.27% (lcr), 90% and 99% 
confidence intervals on the reconstructed (613,6) parameters. The different bins are 
treated as independent. A maximum likelihood fit is performed to the simulated data 
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to extract the signal content, through a global scan in the plane (613, S). The internal 
signal grid used for the fit has a resolution of A813 — 0.1° and A5 = 1°. 

7.3 Combining an OPERA-like detector at 732 km and an 
iron detector at 3000 km 

The numbers of expected events summarized in Tab. 03 and Tab. HI show that the 
sensitivity of the OPERA-like detector degrades considerably at small values of both 
#13 and 5. However, its contribution to the overall measurement is very relevant due 
to the complementary oscillation pattern described in Section |21 Therefore, the effect 
of combining both "golden" and "silver" channels can be already seen for #13 > 1°. 



Confidence level 


"mm 


^max 


68.27 % 
90.00 % 
99.00 % 


60.8° 
49.9° 
29.8° 


98.6° 
106.6° 
118.1° 



Table 7: Allowed regions on the 5 parameter extracted from the analysis of the simulated 
data for 9 13 = 1°,5 = 90°. Both a 5 Eton OPERA-like detector at L = 732 km and 
a 40 Kton iron detector at L = 3000 km are considered. The best fit corresponds to 
13 = 0.9°, 5 = 80°. 

Fig. ITT1 shows the extracted 68.27%, 90% and 99% confidence regions in the (0 13 , 8) 
plane corresponding to four simulated points, if the OPERA-like detector is placed 
at L — 732 km from the neutrino source. Tab. summarizes the projection of the 
confidence level contours onto the 5 axis for 0x3 = 1°, ^ = 90°. It can be noticed that 
all the displayed curves are connected. 

Since the expected background is small (< 10 events) in all bins, the systematic 
uncertainties on this component have a minor effect on the overall sensitivity. The 
same consideration applies to the "silver" signal in the OPERA-like detector. How- 
ever, systematic uncertainties greater than 5% on the golden signal can reduce the 
sensitivity of the iron detector. This is particularly important at small values of 5. 
In our calculation we assumed a conservative 10% uncertainty, which could be further 
reduced by an accurate work on the detector. Tab. IHlshows the effect of this parameter 
on the 99% confidence regions resulting from the simulated point #13 = 5°, 5 = 0°. 

From Tab. El it can be seen that the sensitivity of the OPERA-like detector at 
L = 732 km is limited by backgrounds, mainly from anti-neutrino charm-production 
where the primary lepton is not identified (Section 15.2.2)) . The evaluation of this 
component is based upon the present understanding of the OPERA detector. As 
discussed in Section IPl a detailed dE/dx analysis close to the end point of the tracks 
could further improve the primary muon identificatin by about a factor two. The 
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Figure 11: Expected 68.27%, 90% and 99% confidence regions in the (9is,5) plane 
corresponding to four simulated points. Both an OPERA-like detector at 732 km and 
an iron detector at 3000 km are combined. The stars denote the best fit points, while 
the open circles are the true simulated points. 

corresponding efficiency loss is limited to few percents. When a sizeable "silver" signal 
is expected, such a background reduction would result in narrower confidence intervals, 
as can be seen from Tab. El 
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Table 8: Effect of the systematic uncertainty on the signal in the iron detector for the 
99% CL 5 intervals obtained for 8x3 = 5°, 5 = 0° , when combining the 5 Kton ECC 
detector at L = 732 km and the 40 Kton MID at L = 3000. The best fit for the default 
configuration (i.e., with a 10% systematic uncertainty on the golden muon signal at the 
MID) corresponds to 9 13 = 5.1°, 5 = -2°. 
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68.27 % 
90.00 % 
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82.5° 
65.4° 
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102.6° 
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Table 9: Effect of a background reduction by a factor of two for the OPERA-like detector 
at 732 km and the simulated point # 13 = 5°,5 = 90°. The best fit for the standard 
configuration corresponds to 9x3 = 5.0°, 5 = 86°. 

7.4 Combining an OPERA-like detector and an iron detector 
both at 3000 km 

A location of the OPERA-like detector at 3000 km would reduce considerably most 
background contributions (proportional to 1/L 2 ), with the exclusion of the muonic 
decay of r + events from — > v T oscillations. Tab. El shows the corresponding increase 
in the signal/background ratio for the detection of the "silver" signal events. 

In Fig. [T2l we show the extracted 68.27%, 90% and 99% confidence regions in the 
(9i3,5) plane corresponding to the same four simulated points of Section 1731 In this 
case, the curves extracted for the point 9x3 = 1°,5 — 90° are not fully connected, 
mainly due to the tiny signal contribution expected from the "silver" channel (Tab. EI). 
However, as explained in the previous section, values of 9xs ~ 1° are close to the intrinsic 
limit of the experimental sensitivity, producing wide fluctuations on the number of 
observed events. Actually, by increasing the number of simulated experiments for 
^13 = 1°,5 — 90°, in few cases we observe the presence of disconnected curves also 
with the 732 km baseline for the OPERA-like detector, regardeless of the additional 
background reduction mentioned in Section 17.31 For larger values of #13 the location 
of the OPERA-like detector at 3000 km can produce some reduction in the size of the 
confidence intervals. For comparison, in Fig. ESwe give the corresponding results when 
only the magnetized iron detector at L — 3000 km is used (still considering two beam 
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Figure 12: Expected 68.27%, 90% and 99% confidence regions in the (0 13 ,5) plane 
corresponding to four simulated points. Both an OPERA-like detector at 3000 km and 
an iron detector at 3000 km are combined. The stars denote the best fit points, while 
the open circles are the true simulated points. 

polarities). 

It is worth noting that the scanning load reduces as 1/L 2 and we do expect to scan 
few thousand events in this case. This would allow some potential improvements: 

• remove more than one brick per event. This would allow an increase of the signal 
detection efficiency by about 20% (mainly due to an increase in the brick finding 
efficiency); 

• inclusion of the z7 e — > v T channel when the beam is run with opposite polarity. 
The background is larger by about a factor of two with respect to the u e — > v T 
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Figure 13: Expected 68.27%, 90% and 99% confidence regions in the (0 13 ,5) plane 
corresponding to four simulated points. Only an iron detector at 3000 km is considered, 
with two beam polarities. The stars denote the best fit points, while the open circles are 
the true simulated points. 



channel, while the signal is reduced by the same amount due to the cross-section 
of anti-neutrinos. Since the resulting background will be still small at 3000 km, 
this channel would slightly increase the sensitivity of the measurement in spite 
of the worse signal/background ratio. 

Alternatively, one could also foresee an increase of the mass of the detector and/or 
of the exposure. 
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8 Conclusions 



It was previously shown [23 that looking for "golden" wrong-sign muons at the Neu- 
trino Factory is the most sensitive method to measure simultaneously # 13 and 5. One 
of the main problems concealed in this measurement was, however, pointed out in [2E| : 
due to correlations between # 13 and S in eq. (J2J): degenerate regions in the reconstructed 
(#13, 8) parameter space occur in many cases, thus producing a scenario pretty much 
similar to that of the solar neutrino puzzle before SNO and KamLand. In particu- 
lar, the appearance of different allowed regions in the parameter space, often widely 
separated in the 5 axis, severely reduces the Neutrino Factory sensitivity to the CP- 
violating phase. In [T] it was proposed the use of v e — > v T oscillation to solve the (6* 13 , 8) 
"intrinsic ambiguity" problem, taking advantage of the complementarity between the 
v e v ii and v e — > v T oscillation probabilities (Section EJ). 

Whereas the "golden" channel can be studied using a coarsely grained magnetized 
iron detector [3S] , to profit of the "silver" channel we must use a detector of a different 
kind, capable of separating the "silver" muon signal from the "golden" muon one, by 
means of the different energy distribution of the final muons or by looking for the r 
decay vertex. In P, it was decided to consider the second possibility looking for the 
"silver" muon signal with an OPERA-like detector, making use of the available infor- 
mation on the ECC detector HU] ■ However, since the "silver" channel is strongly 
suppressed with respect to the "golden" channel by the v T N cross-section and by the 
t — ► fi branching ratio (we are indeed dealing with tens of events, at most), extreme 
care must be devoted to the background treatment. This is why we present in this 
paper a dedicated analysis of backgrounds and efficiencies at an ECC OPERA-like 
detector when dealing with the "silver" and "golden" channel at a Neutrino Factory. 
The outcome of this study, presented in full detail in Section is that the three dom- 
inant sources of background to the "silver" channel after emulsion scanning, are (in 
order of importance): wrong-sign muons coming from the decay of charged charmed 
mesons produced in combination with a non-observed right-sign muon; punch-through 
mesons that mimick a charged energetic lepton particle identified as a wrong-sign muon; 
right-sign muons coming from the decay of a r produced through the leading oscilla- 
tion — > u T , whose charge is wrongly identified. Notice that the dominant sources of 
background arise from non-oscillated v e and that undergo CC interactions (charmed 
mesons production) or NC interactions (punch-through mesons): both of them, there- 
fore, decreases like the flux as 1/L 2 and can thus be made negligible by locating the 
ECC detector at a larger distance (as it was indeed the case for the golden signal at 
the MID |2Z|). 

In this paper we show that the intrinsic ambiguity problem is solved for 13 > 1° 
when using a 40 Kton magnetized iron detector, to deal with the "golden" muon signal, 
and a 5 Kton ECC OPERA-like detector, to measure the "silver" muon signal. We 
also include systematic effects in the treatment of the "golden" muon sample at the 
MID, previously not considered. We present (Section^) a refined statistical analysis of 
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the simulated data for a MID at L — 3000 km and an ECC either at L = 732 km or at 
L = 3000 km, taking advantage of the different energy distribution of "silver" signal and 
backgrounds. Below 13 = 1°, the "silver" muon sample at the ECC detector becomes 
statistically negligible. With both the MID and the ECC located at L = 3000 km we 
expect a significant decrease in the background of the "silver" channel, as previously 
mentioned. This translates, for 13 > 1°, in a reduction of the confidence intervals with 
respect to the configuration with the ECC at L — 732 km. 

We have also carried out a rather detailed analysis of the foreseeable scanning load 
at the 5 Kton ECC detector with the considered Neutrino Factory beam (1 x 10 20 useful 
fi + decay in the storage ring), showing that the total scanning load to deal with the 
"silver" and "golden" muon samples and the related backgrounds (~4x 10 4 bricks to 
be scanned in five years) is technologically affordable and does not represent a severe 
problem. 

It must be stressed that in this paper we restricted ourselves to the study of the 
(#13, 5) intrinsic ambiguity, by fixing # 2 3 = 45° and by choosing a given sign for Am^ m 
(in the hypothesis that more information on the three neutrino spectrum will become 
available by the time the Neutrino Factory will be operational). 

We believe that solving the three ambiguities at the same time will need the com- 
bination of different kinds of detectors and baselines: a careful study of the required 
net of detectors is beyond the scope of this paper and it will be presented elsewhere. 
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A Charm-production cross-section 



Neutrino and anti-neutrino induced single charm-production is particularly interest- 
ing to study the strange-quark parton distribution function and the threshold effect 
in the cross-section, associated with the heavy quark production. Over the past 30 
years, many experiments have carried out these studies with complementary tech- 
niques: calorimetry, bubble chambers and nuclear emulsions. These data have been re- 
viewed and combined statistically to extract a world averaged single charm-production 
cross-section induced by neutrino in Ref. [53] • These results have been compared 
with predictions of a leading-order calculation with m c = 1.3 GeV (Fig. (HJ) and the 
agreement is quite good. Unfortunately, there are no inclusive measurements of the 
anti-neutrino induced single charm-production. However, given the good agreement 
between data and theoretical calculations obtained for neutrinos, we use for our calcu- 
lations the parametrization shown in Fig. El which is based on the LO calculation at 
m c = 1.3 GeV of 
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